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Abstract — This letter provides closed-form expressions for the 
outage probability, the average level crossing rate (LCR) and 
the average fade duration (AFD) of a dual diversity selection 
combining (SC) system exposed to the combined influence of the 
cochannel interference (CCI) and the thermal noise (AWGN) in 
Nakagami fading channel. The branch selection is based on the 
desired signal power SC algorithm with all input signals assumed 
to be independent, while the powers of the desired signals in 
all diversity branches are mutually equal but distinct from the 
power of the interference signals. The analytical results reduce 
to known solutions in the cases of an interference-limited system 
in Rayleigh fading and an AWGN-limited system in Nakagami 
fading. The average LCR is determined by an original approach 
that does not require explicit knowledge of the joint PDF of the 
envelope and its time derivative, which also paves the way for 
similar analysis of other diversity systems. 

Index Terms — Level crossing rate, average fade duration, 
cochannel interference, selection combining, Nakagami fading. 



I. Introduction 

DIVERSITY is often used to mitigate fading in wireless 
communication systems and selection combining (SC) 
is one of the simplest diversity combining methods [1]. The 
outage probability (OP) is the first-order statistical measure 
of a diversity system and its typical performance measure. 
The OP of SC systems in an interference-limited environment 
(such as that of cellular mobile systems) for various fading 
channels presented a challenge that was successfully tackled 
by the studies of more then a few researchers, for example, 
for the Rayleigh channel in [2], for the Nakagami channel in 
[3]-[5], and for the correlated Nakagami channels in [6]. The 
OP is also essential to determine other system's performance 
measures, such as the probabilities of bit (symbol) error that 
are also already determined for various modulation techniques 
and corresponding SC receiver structures and fading channels 
[5], [6]. The average level crossing rate (LCR) and the average 
fade duration (AFD) are the system second-order statistical 
measures [7] -[9], which reflect the correlation properties of 
its input/output signals. The expressions for the average LCR 
and the AFD are known for the signal envelope at the output 
of an SC system for independent branch (Rayleigh, Rician 
or Nakagami fading) signals in presence of the thermal noise 
(AWGN) but in absence of the cochannel interference (CCI) 
[10]. In [11], explicit closed-form expressions are derived for 
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the average LCR and the AFD of a dual diversity SC system 
subject to CCI and Rayleigh fading only. In this letter, we 
derive the expressions for the average LCR and the AFD for 
the dual diversity SC system that employs the desired signal 
power algorithm and is exposed to both the AWGN and the 
CCI in generalized Nakagami fading channels. 

II. System and Channel Models 

The dual diversity SC communication receiver is exposed to 
the AWGN and independent and identically distributed (IID) 
single desired signal and multiple CCI signals transmitted over 
Nakagami flat fading channels. The desired signal is received 
over two diversity branches (paths) with same set of average 
fading power fig and fading severity parameter ms, whose 
envelopes xqi~ (1 < k < 2) follow the Nakagami probability 
distribution function (PDF) [12], 

m s \ ms 2x 2ms ~ 1 
fx 0k 0) = ( -^7 ) \ ey P 



r(m s ) 



msx 



x>0, 



where T(-) is the Gamma function, defined by T(a) 



(1) 



l dt [14]. 



The cochannel interference consists of n Nakagami-faded 
IID transmissions (each received over the two diversity 
branches). They are assumed to have same set of average 
fading power £!i and fading severity parameter mi, so the 
envelope of some interference signal i (1 < i < n) received 
in branch k, Wik, is a Nakagami random variable (RV) whose 
PDF is given by [12] 



2w 2mi ~ 
r(mi) 



■ exp 



miw 



,w>0. (2) 



Both ms and mi are assumed to be positive integers. 

In presence of the CCI, the branch selection of the selective 
combiners can rely on three somewhat different decision algo- 
rithms that already received sufficient attention in the literature 
[3, Section IV] and [11]: the desired signal power algorithm, 
the total power algorithm and the signal-to-interference power 
ratio (SIR) algorithm. The desired signal power algorithm 
is the most difficult to implement in practice as it requires 
identification and separation of the desired signal and the CCI. 
However, the known performance analysis of the desired signal 
power algorithm for an interference-limited SC system (the 
outage probability [3] and the LCR and AFD in Rayleigh 
fading [11]) indicates almost identical performance to that of 
the total power algorithm, which is the easiest to implement 
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in practice. For the purpose of analytical tractability, we adopt 
the desired signal power algorithm and reasonably assume that 
the obtained results may also closely indicate the performance 
of the total power algorithm. 

Note that the presence of the AWGN does not significantly 
influence the operation of the desired signal power algorithm 
and its decisions. This assumption is appropriate only if, in 
measuring the largest desired signal power, the average AWGN 
power is taken over a sufficiently long time such that it may 
be considered as a constant across all branches. In this case, 
choosing the signal with the largest desired signal power is 
actually equivalent to choosing the branch with the largest 
signal-to-noise power ratio (SNR). 

Thus, the receiver is assumed to use the SC algorithm 
that selects the branch with the largest short-term desired 
signal power, Xq = max{x^ lf Xq 2 }> or, equivalently, the 
desired branch signal with the largest short-term envelope 
level, xo — maxjxoi, 2:02}- Then, the short-term interference 
power at the output of the SC receiver is consisted of the 
powers of the individual interference signals Wi received in the 
selected branch, and is given by w 2 — w i P> Eqs. (19)- 

(22)] where the PDF of Wi is determined by (2). Considering 
also the presence of the AWGN, the short-term signal-to- 
interference-plus-noise power ratio (SINR) at the output of 
the dual diversity SC system is given by [6] 



max ({sgJjUi) 



y 



(3) 



where a 2 is the AWGN power (being constant over the short- 
term measurement and equal for the two diversity branches) 



and y 



cr is the interference-plus-noise power. For 



an interference-limited system, a 2 = 0, and (3) reduces to 
z\ = Xq/w 2 , denoting the short-term SIR. 

III. Average LCR and AFD 



AFD is defined as the average time that the envelope ratio 
remains below the threshold g, and is defined by 



T g (g) 



F g {g) 
N g (g) 



(6) 



where F g {-) denotes the cumulative distribution function 
(CDF) of the envelope ratio. By introducing (4) into (5) and 
(6), the average LCR and the AFD for the SINR at threshold 
z are found as N z (z) = N g (y/z) and T z {z) = T g (*/z), 
respectively. 



B. Outage Probability 

First, let us determine the PDF of the envelope ratio at the 
SC receiver output. It is calculated from [13, Eq. (2)] applied 
over (4), 



/s(#) = / yfx {gy)f y (y)dy, 

Jo 



(7) 



where f xo (•) is the PDF of the selected desired signal envelope 
xq. It is well known that this PDF can be represented in 
terms of the envelope PDFs and the envelope CDFs of the two 
desired signals as f Xo (x) = F Xol (x)f X02 (x) + F XB2 (x)f Xol (x). 
In Nakagami fading, f Xol (x) — f X02 (x) are determined by 
(1), while the respective CDFs by F Xol (x) — F XQ2 (x) = 
1 — r(mg, msa; 2 /f2s)/r(ms), where r(-, •) is the incomplete 
Gamma function defined by T(a, z) — / z °° t a ~ 1 e~ t dt [14]. 
Thus, the PDF of xq is determined as 



m s 



X 



1 - 



4j.2m s -: 

r(m s ) 

1 



■ exp 



TOgX 



r(ms) 



r m s 



rngx 



2 



x > 0. (8) 



A. Definitions 

To arrive at the desired result, we define the ratio of the 
selected desired signal envelope xq and the interference-plus- 
noise envelope y, 

f~ def Xo 

g = vz = — , 
y 



which will be denoted as the envelope ratio. In absence of 

dcf 

the AWGN, the envelope ratio is represented as g\ = Jz~\ = 
xq/w, denoting the signal-to-interference envelope ratio. 

We will first establish the average LCR of the envelope 
ratio. The average LCR of the envelope ratio at threshold 
g represents the average number of times, per time unit, 
the stationary fading process g(t) crosses threshold g in the 
positive direction, and is mathematically defined by Rice's 
formula [7] 



N g (g) 



gf gg (g, g)dg, 



(5) 



where g denotes the time derivative of g, and f gg (g,g) is the 
joint PDF of RVs g(t) and g(t) in an arbitrary moment t. The 



When ms is a positive integer, T(-, ■) is represented as the 
finite series expansion given by dA.lt in Appendix A. The 
envelope PDF of the output interference signal w is given by 



mi 
o7 



min 2w 2min ~ 



r(min) 



exp 



Ox 



,w>0, 



(9) 

while, by using the standard transformation theory over the RV 
y = \/w 2 + a 2 , the PDF of the output interference-plus-noise 
envelope y is obtained as 



f y (y) = [ g 



x exp 



T{m\n) 
mi . 2 

7^(y 



y > a. (10) 



Applying (8) and (10) over (7), it is possible to obtain the 
closed-form solution for the PDF of the envelope ratio as 
according to the derivation given in Appendix A, which is 



3 



fM 



E 
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r(mg)r(min) 
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(T mi 



1 



(1+ 5 2 //j)™s 



r ! ">s + ^(i + -) 



1 



(1 + 2 5 2/M) ms 



+1 



T 2 



x r(m s +i+j ) ^-(l + 2^-) 



(11) 



where represents the ratio of the average scattered powers 
of the selected desired signal and a single interference signal, 
given by fi = (fismi)/(Qims). 

The PDF of the SINR is obtained directly from (11) 
according f z (z) = f g (y/z)/(2y/z). The CDF of the SINR 
at threshold z is obtained by introducing (11) into F z (z) = 
Jo fg(g)dg, which actually represents the OP of the SC 
system. 

It is possible to determine the closed-form solution of 
the CDF of the SINR, because the first arguments of the 
incomplete Gamma function appearing in (11) are positive 
integers and are expandable as finite series according to ( IA. lb . 
This expansion will yield to simpler definite integrals solvable 
in closed-form involving additional finite series, as follows 



z/n 



a-1 . 

t a -\l + t)- b e- ct dt = {-l) a e c c b - 1 " 

k=0 



(-c) - *[T(* + 1 - b, c(l + z/n)) - T(k + 1 - b, c)] , (12) 

where a and b are positive integers, and c = a 2 m\/Q,\. As 
a result, the closed-form solution of the CDF of the SINR 
involves five-fold summations and will not be presented here 
due to its length. Alternatively, it is possible to numerically 
integrate (11) by using standard numerical techniques and to 
determine the CDF of the SINR for the given set of values for 
z, fi and c. 

For an interference-limited system, (11) specializes into a 
closed-form representing the PDF of the signal-to-interference 
envelope ratio, 



1 



G? 2 /» 



"is -1/2 



y/JIB(ma,min) (1 + j 2 /l') ms+rai " 
1 



x I 



2 + /V.r 



; ms, m s + mjn 



(13) 



which is also derived in Appendix A. In (13), B(-, •) is the 
Beta function B(a, b) = T{a)T{b)/T{a + b), and /(•;•,•) is the 
regularized Beta function defined by I(z;a,b) = I z (a,b) = 
~B z (a,b)/~B(a,b), where B 2 (-, •) is the incomplete Beta func- 
tion defined by B x (a, b) = B(z; a, b) = /„* t a -\l - tf- l dt 
[14]. 

When ms is a positive integer, a closed-form solution for 
the respective CDF is obtainable in terms of the Beta functions 



as given in Appendix A by jA.lj . Consequently, the CDF of 
the SIR at threshold z is given by 



2(-iy 



B(rns,Tnin) \ fi 



B ; ms, 1 — ms — m\n 



B(ms,Tnin) \ 2 



-, x ms — 1 ms — 1 / ■ -, 

1 \ x ^ / TOg + mjn + J — 1 



E 

3=0 



J 



x (^^j B ( _ ^7 ;ms 1 _ m s - min-jj , (14) 

In case of Rayleigh fading (ms = mi = 1), using the general 
identity B(z;l,6) = (1 — (1 — z) b )/b, it is readily verified 
that (14) specializes into a previously obtained result [11, Eq. 
(33)]. 

The results (11), (13) and (14) are new to the best of the 
author's knowledge. 



C. General Expression for Average LCR 

As according to (4), the random process g(t) is related to 
the envelopes of the desired signals 2 fc if) an d the interference 
signals received in the selected branch Wi(t). For an arbitrary 
power spectrum and certain mathematical conditions, each 
envelope's time derivative at any given moment t follows 
the zero-mean Gaussian PDF [7] -[9] and is statistically in- 
dependent from its envelope [8], [9]. Assuming a continu- 
ous wave (CW) transmission and two-dimensional isotropic 
scattering, the Clarke's model determines the variance of this 
Gaussian PDF to be given by <j| = (7r/ m o) 2 ^s/ m s an d 
a ib ~ ( 7r i mi) 2 'Hi /mj for each of the desired signal and the 
interference signals, respectively [9]. It is also assumed that 
the maximum Doppler spreads of the two desired signals 2o/c 
are equal to / m o, and the maximum Doppler spreads of all 
interference signals Wi are equal to / m j. 

To determine the average LCR of the random process g{t) 
by using definition (5), one needs to establish the joint PDF 
of random processes g(t) and g(t) at any given moment t. 
However, we utilize an alternative approach, which circum- 
vents explicit determination of f g g(g,g). From (4), the time 
derivative of the envelope ratio is written as 



1 . x . 1 . g . 

g = - x — oV = ~ x o — y ■ 
y y y y 



(15) 



Conditioned on y, the joint PDF f g g(g,g) is calculated as 



f 99(9,9)= f gg \yig,g\y) f y (y)dy , (16) 
Jo 

where f y (y) is the PDF of the interference-plus-noise envelope 
y given by (10). 

In (16), f gg \y(g, g\y) is the conditional joint PDF of g and 
g given some specified value of the interference-plus-noise 
envelope y, which is expressed as 



fgg\y(9,9\y) = fg\gy(g\9,V) fg\ y (g\v) , 



(17) 



where f g \ y (g\y) is the conditional PDF of g given some 
specified value of y. From (4), f g \ y (g\y) = yf xo {gy), where 
/x (-) is the PDF of the selected desired signal envelope xq. 
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In (17), fg\ gy (g\g, y) is the conditional PDF of g given some 
specified values of the envelopes ratio g and interference-plus- 
noise envelope y. Considering (15), this conditional PDF is 
determined as follows: Conditioned on g and y, g is a linear 
combination of two independent RVs: the time derivative 
of the selected desired signal envelope x (t) and the time 
derivative of the interference-plus-noise envelope y(t). 

We now represent the output interference-plus-noise power 
as a constant plus a sum of m\n IID squared Rayleigh RVs 
Qjj as y 2 = a 2 + Y^i=i a ij> since m i is a positive 

integer. After deriving both sides of the latter expression with 
respect to t, constant a 2 vanishes so that we can apply the 
result [15, Section 3.2.1] and conclude that the time derivative 
of the interference-plus-noise envelope y(t) is a zero-mean 
Gaussian RV with variance equal to the variance of tbi{t) of 
any interferer (a? = cr? .), but only if <t? . is equal for all i, 
as assumed in our case. 

Based on [10, Eqs. (13)-(15)], we also conclude that the SC 
output signal envelope x (t) and its time derivative xo(t) are 
independent RVs at any given t, and that i (t) is a zero-mean 
Gaussian RV with variance equal to the variance of x k(t) of 
any desired signal (cr| o = c| ofc ), since a 2 Qk is assumed equal 
for both k (and x kS are IID RVs at any given t). Consequently, 
g is a zero-mean Gaussian RV with variance 



D. Average LCR for Our System and Channel Models 

Applying (8) and (10) over (21), it is possible to obtain 
the closed-form solution for the average LCR of the envelope 
ratio at threshold g as according to the derivation given in 
Appendix B, and then determine the average LCR of the SINR 
at threshold z as follows 



ms-l/2 



exp(cr 2 mi/ili) 



mjn 

E 

i=l 



m\n — 1 
i-l 



a 2 ra\ 



r(ms)r(min) 
1 



1 a 2 m\ , 

xr|ms+ ^2'^r' ! 



ms-l ^ 



(1 + z/^) m s+*-l/2 

z' 

My 

: 



zjt y 



(l + 2z/^) m s+^i/2 j\ \l + 2z/nJ 



1 cf 2 wli / 2z 



(23) 



For an interference-limited system, (23) specializes into a 
closed-form that represents the average LCR of the SIR at 
threshold z, 



a g\gy ~ 



v 2 a±0 + V 2 



(18) 



Introducing (16) and (17) into (5), and changing the orders of 
integration, we obtain 



N g(g) = gdg fg\ gv (g\g,y) f g \y(g\y) f y (y)dy 

Jo J a 

r-OO />00 

= / fg\ y (g\y) f y (y)dy g f g \ gy (g\g,y)dg . (19) 
Jo Jo 

The inner integral in (19) is calculated by using (18), i.e., 



J g\gy 



POO 

y o gfg\ gy (g\g,y)dg= } 




Substituting (20) into (19) and considering that f g \ y (g\y) = 
yfx (gy), we arrive at the important generalized expression 
for the average LCR of the envelope ratio at threshold g, 



Ng(g) = 



+g 2 <4 r 

2tt J 



fxo(gy) f v (y)dy ■ (21) 



In absence of CCI, Qi = and the PDF of y is given 
by f y (y) = S(y — a), where S(-) is the Dirac delta function. 
In this case, the envelope ratio actually represents the signal- 
to-noise envelope ratio = x /cr, whose average LCR is 
obtained from (21) after applying the fundamental property of 
the Dirac delta function, 



NgM 



fx 0#) • 



(22) 



Introduction of (8) into (22) yields to an expression consistent 
with [10, Eq. (18c)] for the case of dual diversity SC. 



N Zl (z) = V8TrJfl + fi iZ /n 



' T(tos + min — 1/2) 
r(ms)r(min) 



(z/^y 



is-l/2 



' ,W :/„)"■ ^(j-V^ 1 ""-"^ ' I)' 

(24) 

The derivation of the respective LCR of the signal-to- 
interference envelope ratio is provided in Appendix B. For 
the two-dimensional isotropic scattering and CW transmis- 
sions with same maximum Doppler spreads (/ TO j = / TO o), 
a 2 /ct^j. = H and (24) further reduces into 



r(m s + mm - 1/2) (z/M) ms ~ 1/2 

mg+min- 1 



(20) N zi (z) = V8irf m0 



r(ms)r(min) (1 + z/fi) 



X / (2TM7^ ;mS ' mS+min_ ^)- (25) 

In case of Rayleigh fading (m s = m\ = 1), using the 
general identity I(z; 1, b) = 1 — (1 — z) b , it is readily verified 
that (25) specializes into the known result [11, Eq. (32)]. The 
results (23), (24) and (25) are new to the best of the author's 
knowledge. 

Without any loss in generality of the obtained expressions 
for both the OP and the average LCR, it is possible to set 
a 2 = 1 while replacing the average desired signal power 
fis with the average SNR per interferer per branch 73 = 
f2s/c 2 and the average interference power f^i - with the 
average interference-to-noise power ratio (INR) per interferer 
per branch 71 = Qi/a 2 . In this case, fis/^i = 7s/7i an d 
H = (7sTOi)/(titos). 

The AFD of the SINR at the output of the SC receiver at 
threshold z is obtained from T z (z) = F z (z)/N z (z). 
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Fig. 1. Second-order output signal statistics of dual diversity SC system 
with CCI and AWGN vs. SINR threshold for various INRs (71) and fading 
parameters (m = nig = m T ) when n = 3 




Number of Interferers 

(b) Behavior of AFD 

Fig. 2. Second-order output signal statistics of dual diversity SC system 
with CCI vs. number of interfering users for various SIR thresholds when 
mg = mi = 2 



IV. Numeric Results 

Figs. 1 and 2 illustrate the behavior of the second-order 
statistical measures at the output of a dual diversity SC 
system, assuming all CW transmitted signals were exposed to 
two-dimensional isotropic scattering, same maximum Doppler 
spreads (f m i = fmo) an d same fading severity (mg = mi). 

Fig. 1 shows the plots of the normalized LCR N z (z)/f m0 
and the normalized AFD T z (z)f m o in function of the normal- 
ized SINR threshold z/fi (NSINRth) for different values of 
mg and 71 when n = 3. For given ms and 71, the average 
LCR increases with the increase of NSINRth until it reaches 
its maximum at NSINRth= ttiQ, and then decreases (Fig. la). 
Note that th depends on 71 and ms, since n is fixed. As 
expected, the AFD is a monotonically increasing function of 
NSINRth, which is depicted in Fig. lb. 

The average LCR is decreased when the system is exposed 
to a less severe fading channel (i.e. signal envelopes fluctuate 
less rapidly around their means), which is illustrated in Fig. 



la by the increase of ms from 1 to 2. Increasing ms, the 
AFD is only slightly decreased for NSINRth below th , and is 
significantly increased for NSINRth above th as the output 
SINR spends more time around its mean (Fig. lb). 

When NSINRth is fixed to a value below its respective 
tho, the increase of 71 from to 20 dB results in only a 
minor decrease of the average LCR, while the AFD is almost 
unaltered. For values of NSINRth above th , increasing 71 
results in the significant increase of the average LCR (Fig. 
la) and the decrease of the AFD (Fig. lb). 

When the SC system is interference-limited, Fig. 2 displays 
Nzj{z)/ fmO an d Tz I (z)f m Q versus the number of interfering 
users n for different values of the normalized SIR threshold 
z/fj, (NSIRth) when ms = 2. The average LCR can increase 
and/or decrease by varying n, which depends on whether 
NSIRth is selected below or above its th (note that th now 
depends only from n). The AFD is a monotonically increasing 
function of n regardless of the selection of NSIRth (Fig. 2b). 
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V. Concluding Remarks 

In this letter, we provide the analysis of the second-order 
statistical signal properties at the output of the SC systems 
exposed to the combined influence of the CCI and the AWGN, 
which was not previously available in the literature. 

The expression (21) is particularly important since it pro- 
vides solution for the average LCR of the SINR for the SC 
system that uses the desired signal power algorithm in arbitrary 
fading conditions, under the assumption of equal-powered IID 
interference signals in all diversity branches. It incorporates 
the solutions for the average LCRs in the special cases - when 
the SC system is exposed only to the influence of the AWGN 
or only to the influence of the CCI. If the closed-form solution 
of (21) is not obtainable for some specific fading scenario, it is 
possible to apply some known numeric integration technique 
and obtain the desired numeric result with predefined accuracy 
for a given set of input parameters. Furthermore, it can readily 
be proved that (21) is also valid for other diversity techniques, 
such as the maximal ratio combing (MRC) and the equal gain 
combining (EGC). 

Appendix A 

The finite series expansion of the incomplete Gamma func- 
tion •), valid for positive integer values of ms, is given 
by [14] 



r(ms) exp 



Ej_ / m s x 

3=0 J V 



(A.l) 

The finite series expansion of (10) is obtained by using the 
binomial formula, so 

min — 1 / -i \ / \ min 

AW-g C"T )(S) <- 2 » 



2\m\n—i— 1 



x exp 



(m l a 2 \ 2y 2i+1 



■ exp 



my 



\ \l\ J I (min) \ Hi J 

We now introduce (8), (IA. It and (1A.2I ) into (7), change the 
order of summation and integration, solve the inner integrals 
in terms of the incomplete Gamma function r(-, ), perform 
some algebraic manipulations, and obtain 

min— 1 



fM= E 



i=0 



( min — 1 



where 

Mg) 



mi 

n7 



2?y 



2i+2 



■ exp 



r(min) 

4 (g 2 j (m^ min - ' 



miy 
Qi 



(A3) 



fx {gy)dy 



^/JiT(m s )T(min) \fli 



1 



x T I ms + i + 



a 2 mi 



Hi 



(l + <? 2 /V) 



(l+3 2 /M) 



ms+i+l 



i j\ \l + 2g 2 /n) 



(1 + 2g 2 /n)^+^ ~ j 



x T m s + i + j + 1 



a 2 mi 
tti 



(1 + 2g 2 /n) 



.(A.4) 



Expression (11) is obtained directly from dA.3l ) and dA.4| > by 
introducing the substitution for summation index i — > i + 1. 

When the SC system is interference-limited, we apply the 
limes over ( |A.3b and jAAj when a 2 — > 0. In this case, only 
the last term of the sum in (1A.31 >. i — min — 1, remains 
non-zero, while the incomplete Gamma functions in iAAl are 
transformed into the respective Gamma functions, so the PDF 
of the signal-to-interference envelope ratio becomes 
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(A.5) 

Since j is an integer, we can use the identity [14] T(ms + 
min + j) = (mg + min)j r(ms + min), where (-)j is the 
rising factorial (the Pochhammer symbol). We also use the 
finite series expansion of the regularized Beta function /(•;•, •) 
[14] 
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(A.6) 



which is valid when b is a positive integer. Expression (13) is 
obtained by applying dA.61 >. the identity I(z; a,b) = 1 — 1(1 — 
z; b, a) and the definition of the Beta function over ( IA.4b . The 
respective CDF of the signal-to-interference envelope ratio is 
obtained by integrating dA.51 > according F gi (g) = f gi (t)dt. 
Thus, after changing the order of integration and summation, 
performing some algebraic manipulation and applying the 
definitions of the Beta function and the incomplete Beta 
function, we obtain 
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Appendix B 

In order to derive (23), we introduce (8), ( lA.lb and dA.2b 
into (21), change the order of summation and integration, solve 
the inner integrals in terms of the incomplete Gamma function 
r(-, •), perform some algebraic manipulations, and obtain 
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where 



mi 



2?y 



2i+l 



■ exp 



r(min) 

£lg r(mg)r(min) \ f2i 
x T ^m s +i+l/2, 
1 



fx (gy)dy 
l 



(1 -|- 91)m s + i+l/2 



a 2 mi 

ms — l 



(l+3 2 /^) 



(l + 2 5 2 /^) m s+*+i/2 j\\l + 2g 2 /n 

x r (m s + i + j + 1/2, (1 + 2.g 2 / M ) 



(B.2) 



Expression (23) is obtained directly from (IB. Il l and JB.2t by 
introducing the substitution for summation index i — ► % + 1, 
and then setting N z (z) = N g (^/z). 

When the SC system is interference-limited, we apply the 
limes over ( IB. Il l and (IB.2l i when rr 2 — > 0. In this case, 
only the last term of the sum in ( IB.U . i = m\n — 1, 
remains non-zero, while the incomplete Gamma functions in 
(IB. 2b are transformed into the respective Gamma functions, so 
the average LCR of the signal-to-interference envelope ratio 
becomes 
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(B.3) 

Since j is an integer, we can apply the identity involving the 
Pochhammer symbol [14], T(ms + m\n + j — 1/2) = (mg + 
m\n— l/2)j r(mg +min— 1/2). Expression (24) is obtained 
by applying JA.61 , the identity /(z; a, 6) = 1 — 1(1 — z; 6, a) 
over (IB.3K and then setting N zi (z) = N gi (y / ~z). 
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